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LIFT AND 
EDGES 


MCMENT ON THIN ARROWHEAD WINGS WITH SUPERSONIC 
OSCILLATING IN SYMMETRIC FLAPPING AND ROLL 
AND APPLICATION TO THE FLUTTER OF AN 
ALL-MOVABLE CONTROL SURFACE 


By H. J. Cunningham 


SUMMARY 


A theoretical treatment is presented for determining lift and moment 
on thin arrowhead or pointed-tip wings of which the delta plan form is a 
special case with an unswept trailing edge. On the basis of linearized 
supersonic potential -flow theory for the symnetric flapping and rolling 
modes of harmonic oscillation, expressions have been developed for the 
section and total forces and moments to the third power of the reduced 
frequency. A limitation to the Mach number range is that the component 
of flow normal to all edges is supersonic or sonic. 

Semple results are given in the form of curves to show that retention 
of terms to the third power of the frequency gives good accuracy over a 
range of frequency that covers many practical flutter applications. Plots 
are given of the spanwise distribution of section force and moment for 
symmetric flapping and rolling modes . Approximating the section force 
and moment by multiplying the section quantities of a rigid translating 
wing by the flapping and rolling mode shapes (termed a "finite -wing strip 
theory") is shown to result in an overestimation of forces and moments in 
comparison with the results of the present analysis. 

A modal type of flutter analysis for an all-movable control surface 
wherein all bending and twisting flexibilities are effectively concen- 
trated in a supporting shaft is made by use of natural (coupled) modes. 

A sample flutter analysis using two coupled modes is given for a Mach 
number of 1.6 for an arrowhead wing with the leading edge swept back 
45° and the trailing edge swept forward 15°. A near coincidence of the 
first two natural frequencies is found to be detrimental as regards 
stiffness required to prevent flutter. 
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INTRODUCTION 


In the continuing problem of analyzing the flutter of wings and 
control surfaces on high-speed aircraft, the flutter of all-movable con- 
trol surfaces is of current concern. Inasmuch as such all -movable sur- 
faces are controllable in pitch about an axis, their root support may 
not be so rigid as that of a fixed-root airfoil. Consequently, the 
natural vibratory modes and the flutter mode may involve appreciable 
amounts of symmetric flapping and roll of the elastically undeformed 
surface about a relatively flexible root support. In order for a flutter 
analysis to be carried out, the distribution of air forces on the har- 
monically oscillating surface is needed. 

The treatment of the present analysis is limited to arrowhead plan 
forms with supersonic leading and trailing edges. Expressions are pre- 
sented for the velocity potential and the distributions of associated 
forces and moments . The velocity potential is obtained by expanding the 
general velocity potential for "purely supersonic" plan forms, developed 
in reference 1, in terms of a frequency parameter, as was done in refer- 
ence 2. The expansion is carried to the third power of the frequency, 
inasmuch as expansions of this order are considered sufficient for a 
large number of practical applications. In reference 2 the arrowhead 
plan form was assumed to be oscillating in. vertical translation and pitch 
and, in the present study, the plan form is assumed to be oscillating in 
a symmetric flapping mode and in roll. The velocity potentials for 
flapping and roll are used to obtain force and moment coefficients for 
strearnwise sections of the wing. Total force and moment are also 
obtained. 

Several other investigations have dealt with the forces on a rolling 
arrowhead wing with supersonic edges. In reference 3 an expression is 
given for the total rolling -moment coefficient for an oscillating delta 
wing. This rolling moment is exact within the framework of linearized 
supersonic potential -flow theory and is used in the present paper to 
illustrate the accuracy of the total rolling moment approximated by the 
frequency-expansion method for the special case of a delta wing. Refer- 
ence 4 includes expressions for the clynamic stability derivative for 
dam ping in roll of the arrowhead wing based on retention of only the 
first term in the frequency expansion for the velocity potential. 

For the purpose of illustrating a main use of the section forces and 
moments for flapping motions, a coupled-mode type of flutter analysis is 
outlined for an all -movable control surface. Each coupled mode is 
assumed to contain components of flapping, vertical translation, and 
pitching. A numerical example of such a flutter analysis is included. 
Other uses for the section forces and moments for the flapping mode and 
some properties of section and total force and moment coefficients are 
also considered. 
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SYMBOLS 


a 

a i 

A i 

Af, Af *, Aj_f, 

A lcp 

b 

C = cot A 
^mn 



x i 

f, h, <p 


F e , H e 
— mu ’ — mn 


&, S ± 


velocity of sound in main stream, ft/sec 

functions of x, 35, and M defined following equation (4) 

generalized mass parameter, A* Qpl? (See equations (21) 
and ( 27 )) 1 

generalized mass terms defined by equation ( 21 ), slug-ft^ 

Bjjf, Bgf coefficients in velocity potential for flapping 
wing defined by equations (A4) to (A 7 ) 

coefficients in velocity potential for rolling wing defined 
by equations (A9) and (A10) 

semlchord of wing at root or at midspan, ft 


generalized aerodynamic force coefficients defined by 
equations (24) and ( 25 ) 


functions of x-j_, y, and p defined in equations (Bl) 
and reference 2 

angular flapping displacement of wing semispan about axis 
y = 0 , positive for wing tip down 

amplitude of flapping displacement f, radians 
component of flapping displacement f in mode i, radians 

time derivatives of f, h, and cp, respectively 
functions of x^_, y, and p defined by equation (B9) 

structural damping coefficients, general and in mode i 
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h 


vertical displacement of wing at x = xq and y = 0, 
positive downward, ft 

amplitude of vertical displacement h, ft 


component of vertical displacement h in mode i, ft 


I a , section mass moment of inertia in pitch about x = xq 

and in flapping about y = 0, respectively, slug-ft^/ft 

J mn function of Xj_, y, and 0 defined in equation (B2) 


k reduced frequency, boi/V 

Lf, components of section lift and section pitching -moment 

coefficients due to vertical translation and pitching 
of rigid wing, defined in reference 5 (i = 1, 2, 4) 



+ iL 


2,f 


complex section lift coefficient due to symmetric flapping 
oscillation 


L l,cp + iL 2,cp complex section lift coefficient due to rolling oscillation 

Lj_, M ± , f , f total coefficients, spanwiBe integrals of section 
J ;X coefficients under bar 


m 

M 


section mass, slugs/ft 
Mach number, v/a 


M l,f + 

M l,cp + ^2,cp 
M i,f’* M i,cp* 

\l + “^,2 


section and total pitching moment about Xq, positive 
leading edge up, ft-lb/ft and ft-lb, respectively 

complex section pitching -moment coefficient about Xq due 
to symmetric flapping oscillation 

complex section pitching -moment coefficient about xq due 
to rolling oscillation 

parts of and obtained when Xq = 0 

total rolling -moment coefficient due to rolling oscillation 


Ap 


local perturbation pressure difference, positive up, lb/sq ft 
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P, P section and total lift force, positive down, lb/ft and 

lb, respectively 

P _ section cross product of inertia referred to x = Xq, 

y = 0, slug-ft 2 /ft 

Qj_ generalized aerodynamic force in mode i, Ib-ft 

R = fir- j ) 2 - p 2 (y - n ) 2 

s *> s f section static mass unbalance about Xq (positive leading 

edge up) and about y = 0 (positive right tip down), 
respectively, slug-ft/ft 

function of X]_, y, and p, defined in equation (B2) 


t 


time, sec 


u, Uq, v, Vq characteristic coordinates defined in appendix A 
V velocity of main stream, ft/sec 


w 


5W 


downwash at surface of wing; that is, increment of 
vertical velocity imparted to air stream by position 
or motion of wing surface, positive when air is deflected 
away from that surface, ft/sec 

virtual work, see equation (22), lb-ft 


x', y', z' rectangular Cartesian coordinates, ft 

x, y rectangular Cartesian coordinates, nondimensional in terms 

of chord 2b: x = — y = 

' 2b 2b 


*0 


X 1 

y t 

Zm(x,y,t) 


x-coordinate at pitch axis and coincident pitching-moment 
axis 

x-coordinate of trailing edge 
y-coordinate at wing tip, C/(l - CD) 

vertical displacement of wing at point (x,y), positive up, 
ft 
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a 

a 0 


angular pitching displacement of wing about Xq, positive 
leading edge up, radians 

amplitude of pitching displacement ct, radians 



_ 1 - pc 
1 + 3C 

Ih half -span integrals of lift and moment coefficients, 

defined by equations (2 6) 

9 ^ component of pitching displacement a in mode i, radians 


^(t) 


ii 


sweep angles of leading edge and trailing edge, respectively, 
positive for sweepback, deg 

rectangular Cartesian coordinates used to represent location 
of sources in xy -plane 

generalized coordinate indicating amount of mode i con- 
tained in a given displacement 

amplitude of generalized coordinate 
air density in main stream, slugs/cu ft 


= p2c 2 - 1 


9 

^0 

0 

OJ 

55 


rolling displacement of wing about y = 0, positive right 
wing down, radians 

amplitude of rolling displacement 9, radians 

disturbance velocity potential, ft^/sec 
circular frequency of oscillation, radians/sec 
circular natural frequency of mode i, radians/sec 
frequency parameter, 2kM 2 /p 2 
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j j absolute value of quantity between vertical bars 

A bar over a quantity indicates value of quantity integrated over 
the span (except 05) 


AMLYSIS 


Equations for Velocity Potentials 

Integral form .- On the basis of a rectangular coordinate system 
moving forward with the wing at a constant supersonic speed V in the 
negative x' -direction (see fig. 1), the differential equation for the 
propagation of small disturbances to be satisfied by the velocity 
potential 0 is 


U 8 

a 2 V3t 


+ V 



I— ¥— + — — aU u 

\ 8 ( x ') 2 8 ( y ') 2 8 ( z ') 2 / 


( 1 ) 


The principal boundary condition that the velocity potential must 
satisfy is that the flow at the wing be tangent to the wing surface. 
This boundary condition can be expressed by 



where is the vertical displacement from z' = 0 of any point of 

the wing (fig. 2) and w is the downwash; that is, w is the vertical 
velocity imparted to the air stream by the position or motion of the 
wing surface, positive when directed away from that surface. In accord- 
ance with linearized theory this boundary condition is evaluated at the 
projection of the wing on the z' = 0 plane, and thus no effect of 
wing thickness is taken into account. 

The solution for 0 that satisfies the present boundary-value 
problem (wing with all supersonic edges ) is s imi lar to the solution 
developed in reference 1 for an infinite swept wing. The problem is 
satisfied by a distribution of sources on the upper (z = +0) wing sur- 
face and, in the absence of thickness, the same distribution but with 
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opposite sign on the lower (z = -0) surface. At this point in the 
development of 0 and for the rest of this analysis, it is convenient 
to use the coordinates x,y as nondimensional quantities, obtained by 
dividing x',y' by the root chord 2b. Then, in a form slightly dif- 
ferent from that appearing in equation ( 9 ) of reference 2, 0 can be 

expressed as 


0(x,y,+O,t) = 2k )% os f5£l\ 3 ^ ( 3 ) 

Jt R \ M / 


where 


s = » 
P 2 


(with k = fep), R = - t) 2 - P 2 (y - T]) 2 , and g 


and t], which denote the location of sources, are nondimensional vari- 
ables of integration in the x and y directions, respectively. The 
region of integration S is, for example, the shaded area in figure 5(a) 
prescribed by the forward facing Mach cone with apex at the point 
(x,y,+0) . The local downwash w is positive when the air is diverted 
outward from the surface by the inclination or motion of the wing. 


If the method employed in reference 2 is followed, the integrand of 
equation (3) is expanded in a Maclaurin's series with respect to m. 

When the series is limited to the third power of and the terms are 
collected with respect to | , the result is 


where 


0(x,y,+O,t) 



S 

a 5 ^ + bQR + b^R^dri dg 


a 0 



+ i 


x 3 

6 


a l 


8. g 


a 5 


ici) + a^x - i ^x 2 



-i 


*3 

T 


00 
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■u _ a 2 , < a? 

Do 2 + 1 ~ 

ar 2 m* 

*1 = -1 


2M 2 


A difference in the present notation from that of reference 2 is that 
herein % , r\, and 35 and consequently R, a n , and b n are used in 

nondimens ional form. 


Reference 2 treated the vertical translation and pitching modes of 
oscillation. The present paper treats the symmetric flapping and rolling 
modes illustrated in figure 3, and the integrated expressions for 0 for 
these modes are given in the following sections . 


Integrated form for symmetric flapping and roll .- The displacement 
and the resulting downwash w on the upper surface in equation (2) are: 


For f , 

= -2bf |y |' 

r 

v = -2bf jy| 

For cp, 

^ = -2bcpyl 
v = -2b$y I 


(5a) 


(5b) 


Evaluation of the integrals indicated in equation (4), after sub- 
stitution of the appropriate downwash, results in the following forms 
for the velocity potential on the upper surface: 




- p^y 2 + A^* tanh“ B 


- ikl + 

X + p|y | 


B lf 


cos-1 

p(Cx + y) 


+ Bgf 


cos 


-1 


x - p 2 Cy 
p(Cx - y) 


(6a) 
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0. = l£(2b) 2 


A m Jx 2 - (J^y 2 - Btj. coS -1 x + + 

‘Pv K 17 IT . p ( Cx + y) 


B 2 , cos' 1 V P - c r T 
2f p(Cx - y) 


(6b) 


where, with reference to the shaded regions of integration of figure 3* 
the terms involving the coefficients Af and A^ contain contributions 

from regions I, II, and III, and the terms with the coefficients A^*, 
B-j_£, and BQf are contributed solely from regions I, II, and III, 
respectively. The coefficients Af, A^*, A^, B^, and Bgj are given 

in appendix A. Appendix A also indicates how the integrations of equa- 
tion (4) were carried out, as in reference 2,vby the convenient use of 
the characteristic coordinate Bystem shown in figure 3(b) . The tanh - ^- 
function with its coefficient Af* arises as a result of the disconti- 
nuity of spanwise slope at y = 0 in the symmetric flapping mode. 

For the limiting condition of sonic leading edges with pc = 1 
(that is, with the Mach lines from the leading apex coincident with the 
leading edges and the entire plan form contained in region I of fig- 
ure 3), eq.ua tions (6) reduce to the forms 

(*4) S C=1 - 1r (2b)2 (W* 2 - U 2 * 2 + V <Ta) 


(Hoi = f< 2b)2A i <eF - (7b > 

The coefficients A lf and A 1( p and the development of equations (7) 
are also given in appendix A. 

Concerning the development of velocity potentials, it is to be 
noted that thus far the only restriction with regard to the trailing 
edge is that the component of flow normal (in the xy -plane) to all 
points of the trailing edge is supersonic or at least sonic. Thus, 
the trailing edge can be jagged or curved - or need not even be sym- 
metric about the midspan. 
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Section Forces and Moments 

The foregoing expressions for velocity potential are now used to 
obtain the forces and moments on any wing section y, as in figure 1 , 
and the observation is made that only those restrictions on the 
trailing edge discussed in the preceding paragraph for the velocity 
potential apply also to the section forces and moments. 

The difference between local pressures on the upper and lower 
surfaces of the wing is 




(8) 


where p is the density of the main stream and 0 is the velocity 
potential for the tipper surface, given in equations ( 6 ) or ( 7 ). With 
the signs used, tip is positive upward. 


The expression for 
figure 1 is obtained by 


the force, positive downward, on section 
integrating as follows: 


y 


of 


P = 



hp dx 


(9) 


where xq(y) is the x coordinate of the trailing edge. The section 
pitching moment (positive leading edge up) about an axis at x = Xq is 




Jy/c 


(x - xq)Ap dx 


( 10 ) 


Substitution of equation ( 8 ) into equations (9) and (10) yields 


P = 2pV 


y/c 


+ 2ik 


X 1 

,y/c 


0 dxi 


( 11 ) 


Ma = MV 


x0 


xi nx x 

y/c Jy/c 


0 dx + 2ik 


/y/c 


c0 dxi 


- 2xobP (12) 
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Upon substitution of 0 from equations (6) or (7) into equa- 
tions (ll) and (12), the section lift and pitching moment can be 
written as 

P = -4pbV 2 k 2 ei^t [f 0 ( Lljf + ±L 2}f ) + cp 0 (lrj_ j( p + iL^)] (13) 

% = -4pb 2 V 2 k 2 e^[f 0 (M ljf + iM 2>f ) + cp 0 ( M l,cp + ^ 2 ,^ W 

In the past, some writers (for example, ref. 2) have used the nota- 
tions In, L 2 , M]_, and M 2 as flutter coefficients associated with 

bending or vertical translation of wings . Inasmuch as each section of 
the wing translates vertically in the two modes of oscillation treated 
in the present paper, the additional subscripts f and cp are used in 
equations (13) and (14) for the purpose of differentiating these new 
coefficients from those for a two-dimensional wing or a wing translating 
vertically as a rigid unit. 

The section contribution to the bending moment about the wing root 
as a result of symmetric flapping (positive in the wing -tip -down direction) 
is 


2b|y|p = -4pb 2 V 2 k 2 e icat f 0 (2|y|L 1>f + i2jy|L^ f ) (15) 

and the section contribution to the rolling moment due to roll (positive 
right wing down) is 

2byP = -4pb 2 V 2 k 2 e i “ t cp 0 (2yL 1 ^ q) + iPyLg^) (l6) 


The section force and moment coefficients and of equa- 

tions (13) to (l6) are given in appendix B. 


Total Forces and Moments 


The total aerodynamic forces and moments needed for a flutter 
analysis or for any other use can be obtained by integrating the section 
quantities over the span. At this point in the present analysis further 
treatment is restricted to plan forms with a symmetric trailing edge and 
particularly to a swept trailing edge which meets the leading edge at a 
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pointed tip as illustrated in figure 1. The angle of sweep of the 
trailing edge is Amg (sweepback is positive), and at the positive 

wing tip the coordinate y = y^ = w here D = tan 

Integrating equations (13) and (14) across the span yields 


P = / 2bP 
J Span 


dy = 2 2bP dy = -8pb 2 V 2 k 2 e iajt f f + iL 2 ^ f ) (17) 


Ma = / 2tM a ^ = 2 j 1 20^ dy = -Qpb^V^e^fJ^ f + lMg f ) (l8) 

JSpan JO X ’ 


where the contributions due to roll are zero by virtue of their anti- 
symmetry. The total root bending moment, obtained by integrating equa- 
tion (15) over the positive half -span, is 


ny t P y t 

^b 2 |y|P dy = -8pb3v2k2ei^tf 0 PyfLn f + iL 2 f W (19) 

Jo Jo K J 


and the total rolling moment, the full -span integral of equation (16), is 


2 jT Yt lb 2 yP dy = -8pb5v 2 k 2 e imt q) 0 2 JJ* + ±L 2 , cp )dy 


" -SpbW^^M^ + iM 9j2 )dy 


( 20 ) 


Bars over the coefficients in equations (17) , (18), and (20) denote 
total coefficients; for example 

L i,f = 2 hyf dy 

Expressions for total force and moment coefficients are given in 
appendix B. 
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Analysis of Flutter Involving Flapping and Pitching 

On the basis of the simplifying assumption that all flexibility is 
concentrated in the supporting structure at the root, a method is out- 
lined in this section for analyzing the flutter of an all -movable con- 
trol surface . The support flexibility allows the control surface to 
pitch, to flap, and to translate vertically". Coupled modes are used 
because experience has shown that the natural vibration modes of the 
control surface are often strongly coupled, and the necessity for 
obtaining frequencies of hypothetical uncoupled modeB is thereby 
eliminated. 

The present analysis is adapted from that of reference 6 and only 
the essential points are given herein. The use of two coupled modes is 
illustrated, but any number can be employed. For two inodes 

h(t) = H 1 6 1 (t) + h £ g 2 (t) 
a(t) = 0-j_|-^(t) + Sgl^t) 

f(t) = f]_£]_(t) + f 262 ^^ 


where h, a, and f are the displacements in vertical translation, 
pitching, and flapping, respectively j hj_, 0 and f.j_ are the components 

of the three types of motion in mode i and are not functions of the 
spanwise coordinate in the present usage; £^(t) = is the gener- 
alized coordinate in mode i and is its amplitude, which In general 

is complex to account for phase differences of the modes. 


The generalized mass of the semispan for mode i is 


A ± 



Iffj.* + 2 


( s a ] 


h i e i + 






dy> 


( 21 ) 


where 


m 

wing mass per 


pitching mass 

If 

flapping mass 


unit Bp an 

moment of inertia about 
moment of inertia about 


x = Xq per unit span 
y = 0 


per unit span 
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S a static mass unbalance about x = xq per uni t span 

Sf static mass unbalance abbut y = 0 per unit span 

p af mass product of inertia about x = x Q , y = 0 per unit sp an 

The virtual work done on the semispan control surface as it moves 
through the (real) virtual displacements Bg and 5S 2 is 


sw = 511 fV + Ma01 + pfi f y, i) dyi 

S *2 + M a 9 2 + »2|y’l) 


= + ^2^2 


( 22 ) 


where 

P = -UpbV^e 1 ^ 


Zo, 

b 


( L i + iL 2 ) + + ^4) + ^(^f + iL 2 ,f) 

M a = -^pb 2 V 2 k 2 e 1 “ t |^(M 1 + iM 2 ) + cx 0 (m 5 + IM^) + f 0 (M 1#f + iM^ f ) 


(25) 


^ 8ft ± 

The generalized forces , after expansion and collection of terms, can 
be represented as 


Ql - Spb^tu^e^fc (|_i C i + L2C12) 
0,2 — 8pb5 l m 2 e^ a) ^^|j_C2]_ + |_p Cpg^ 


( 24 ) 
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in which, by matrix notation, 


r 

-Cu 


U) \ b b f i b e i e i f i 

_C 12 

< 

► = 

fl h 2 h,e 2 fi0£ v e Va .A 0i9a Va 

“ C 21 


"t 1 “t 1 '#1 T 1 *A 9 f 1 9 2 9 i 

-c 22 


( h 2\ 2 f 2 h 2 h 2 0 2 .p q h 2 f 2 ^ 2 0 2 h 2 a 2 n ^ 

U ) b b f2 ® 2 b f2 b 92 9 2 f2 

- 


A 1 

r 2 


r 4 

<r (25) 
r 6 
r 7 

r 8 

l r 9j 


The minus signs appear with the generalized aerodynamic-force coeffi- 
cients C mn 80 that their sign in the flutter determinant to follow 
coincides with that of reference 6. As used herein the elements C. 
are in a nondimens ional form, however, and differ in that respect 
from reference 6. 


mn 


The integrals are 


py+ 


r l " 


(l x + iL 2 )dy 


r 3 ■ JjXb + 


5 t ( Ll ^ f + ^f) 6 ^ 


r = = 


r 4 = 


r 6 " 


r ? = + w 2 )<3y 

r 9 = f t ( M l,f + iM 2,f) <3 ^ 


= 2 |y|( L l + i 4 2 )dy 

= Jj* 2 |y |( l 3 + iL 4) 

= 2 |y |( L l,f + iL2,f)<iy|-(26) 

= + ^4)^ 
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where y + = — . Expressions for the section coefficients L* and 

z 1 - CD 1 

Mi (i = 1 , 2 , 3, k) are given in reference 2 and their spanwise 
integrals Lj_ and Mj_, in reference 5. It is to be noted that the 
integrals I\j_ are, for convenience, evaluated over the half -span; 
whereas the barred quantities Li and Mi represent the full-span 

integrals. Accordingly, one -half of the latter quantities are to be 
used. The integrals and equal one -half the full -span or 

total coefficients appearing in equations (17) and (18), respectively. 
The integral Tg appears on the right side of equation (19) and its 

integrand is given by equation (BIO) of appendix B. The integrands of 
T2 and are obtained by following a procedure parallel to that 

described for obtaining the integrand of Tg beginning with 
equation (b8). 


Finally, with the intermediate steps of reference 6 omitted, the 
equations representing the dynamic equilibrium of flutter axe 


- (^) 2 (l + i gl ) + 

6piWetat + jU . (*)*( X + ig£ ) 


+ ^ 2 C 12 “ 0 


+ C 


221 


= 0 


(27) 


where 


a i 


Aj 

8pb5 


In equation (27) a^ and C TOT1 are nondimens lonal for convenience and 

only in this matter of form and the number of modes treated does equa- 
tion (27) differ from equation (12.30) of reference 6. 


The nontrivial flutter solution is obtained by requiring the 
determinant of the matrix of coefficients of ^ to vanish as follows 

2 


a n 


&) i 1 + lg i) 

C 21 


+ C 


11 


'12 


a 2 


1 - (irjT 1 + ig 2) + °22 


= 0 


( 28 ) 


A sample flutter solution is carried out in- appendix C. 
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DISCUSSION 


Some Characteristics of Forces and Moments 

Spanwise distribution of forces and moments .- The distribution over 
the span of the various components of lift and moment is examined as a 
matter of interest. As an example the coefficients of equations (15), 
(l4), and (15) have been calculated at a number of spanwise stations for 
a 45° delta plan form (C = 1, D = 0) with Xq = 0.5, M = 1.6, and 

k = 0.5 • These results are plotted as functions of the spanwise coordi- 
nate in figure 4. For the symmetric flapping mode the section lift and 
pitching moment are shown in figures 4(a) and 4(b), respectively. For a 
rolling oscillation the section lift is shewn in figure 4(c). In each 
of these figures the solid-line curves result from the inclusion of all 
terms to the third power of the reduced frequency k, and the short-dash 
curves result from the inclusion of terms to the first power of k. As 
can be noted from the proximity of the solid-line and short-dash curves, 
the addition of the second term with the next higher power of k to the 
first term in each coefficient has a moderate rather than a large effect 
everywhere on the span for this sample set of parameters, even for the 
rather high value of k = 0.5« 

In flutter analysis a fairly common practice, when the exact span- 
wise distributions of force and moment are not known, is to use the 
spanwise distributions for rigid-wing motions multiplied by the mode 
shape (a procedure referred to as strip theory) . This method was used, 
for example, in reference 9 In analyzing the flutter of cantilever wings 
of rectangular plan form. Therefore, as an interesting comparison, the 
results of a "finite -wing strip theory" are also shown in figure 4. 

These curves were obtained by multiplying the section coefficients (to 
the third power of k) for a rigid, vertically translating wing (with 
the same plan form and stream flow parameters) by the appropriate 
deflection mode shape, 2ly| for symmetric flapping and 2y for 
rolling. The factor of 2 appears with the strip-theory results because 
the unit of vertical, translation is the root semichord b, whereas the 
flapping and rolling deflections axe referred indirectly to the full 
chord 2b. The strip -theory forces and moments are moderately to 
appreciably different from the results of the present method, particu- 
larly in the case of the in-phase components L^f, M-^^, and Li,q>- 

This difference is a direct measure of the effect of aerodynamic induction 
associated with the flapping and rolling mode shapes treated in the pres- 
ent paper. 

Accuracy as function of reduced frequency .- The coefficients of the 
present analysis have been obtained from a power-series expansion in 
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terms of the reduced frequency k. An evaluation of the loss of accu- 
racy as k increases from zero is necessary to establish confidence in 
the use of the coefficients and to ■ determine an upper boundary of k 
beyond which the present coefficients, which extend to the third power, 
should not be used. 

As an example, for a Mach number of 1.6 and with Xq = 0.5> the 
total lift coefficients L]_ f and L2 f of equation (B27), the pitching- 
moment coefficients M^_ ^ and °f equation (B28), and the rolling- 

moment coefficients ]_ and 2 of equation (B29) have been evalu- 
ated and the expressions for delta wings with supersonic and sonic edges 
are as follows: 


= 0.171077 - 0.154852k 2 + 0.0568936k 14 - - 0.0155879k 6 


= 0.533761^ - 0. 168442k + 0.0953473k 5 - 0.0312471k 5 


M- 


= 0.102646 - 0.096323k 2 + 0.044250^ - 0.0127538k 6 


M. 


= 0. 266880^ - 0; 112297k + 0.071511k 3 - 0.0249977k 5 


= 0.136862 - 0.0770583k 2 + 0.0252860^ - 0.00566834k 6 

c 5 


“9,2 


= 0.533761^ - 0.112297k + 0.0476736k 5 - 0.0124988k 5 


(29) 


As explained in appendix B the total coefficients for the special case 
of the delta plan form are obtained to the seventh power of the frequency 
by application of the reverse-flow theorem. 

The coefficients of equation (29) are plotted in figure 5 in such a 
way as to indicate the use of only the first term (the horizontal line 
in' each case), then the successive addition of the second, third, and 
fourth terms in each coefficient in order to demonstrate the corre- 
sponding progressive improvement in accuracy of these series for higher 
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frequencies. The highest power of k contained in each calculation is 
indicated near the right end of the curve. For the rolling -moment coef- 
ficients of figure 5(c), the exact (not frequency expanded) result of 
Miles (ref. 3 ) for delta wings is also shown and the convergence of the 
series result is evident. (The exact result was obtained by using in 
equation (70b) of reference 3 the aerodynamic functions fx tabulated 
in reference 7.) The relation between the present notation and that of 
reference 3 is 


2C + uya) 

2 P c3 


(30) 


where C-j is the 

‘■p 

respect to rolling 
tive C 7 of, for 

of Mq ^2 with the 


derivative of rolling -moment coefficient with 

velocity (dynamic -stability notation). The deriva- 
example, reference 4 corresponds to the first term 


relation 

-kJL 2 (1 - CD)3 

2C 7 = — 3U£ 

Z P c 3 


(31) 


If an inaccuracy of 10 percent can be tolerated, the first two 
terms of equations ( 29 ) are acceptable over a range of 55 frcan 0 to 2.2 
for the coefficients of figures 5(a) ana 5(' b ), and from 0 to about 2.7 
for the coefficients of figure 5(c). Over these ranges of 55 the curves 

utilizing as the highest power k^ and k® are essentially exact, and 
in figure 5(c) the exact result itself from reference 3 is available for 
comparison. 

The accuracy of the curves of figure 5 extends to higher values of 
55 than was noted in reference 5 Tor delta wings oscillating in vertical 
translation and pitch, and also in reference 8 for two-dimensional wings. 
The reason for the improved accuracy appears to be that, since the local 
amplitude of the oscillation modes treated herein (symmetric flapping and 
rolling) is zero at the wing center section and increases progressively 
toward the wing tips, the wing section that might be termed the "repre- 
sentative section" is located farther outboard where the local reduced 
frequency is smaller than the basic or reference k, which is defined in 
terms of the root chord, k = o5b/v. 

The accuracy shown by figure 5 applies specifically to delta wings. 
Some information concerning the accuracy for the general arrowhead plan 
form is shown by table III of reference 5 and the associated discussion, 
where sweepforward of the trailing edge (decrease in aspect ratio) is 
shewn to result in improved accuracy, and conversely, sweepback of the 
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trailing edge (increase in aspect ratio ) results in decreased, accuracy 
for the total coefficients and for translation and pitching 

modes of oscillation. In fact when Aijg was decreased from 10° to -20°, 

the accuracy appro xima tely doubled. For the modes of oscillation treated 
in the present paper no reason for expecting a different trend of accu- 
racy as a function of the trailing -edge sweep angle is apparent. 


Results of a Sample Flutter Calculation 


In the analysis a method is outlined for analyzing the flutter of 
an all -movable control surface with the flexibility concentrated in the 
root support so that the motions of the surface are made up of flapping, 
pitching, and vertical translation. A numerical example of a flutter 
analysis for the configuration shown in figure 6 is carried out in 
appendix C, where the pertinent parameters and mode shapes are given. 

The results are plotted in figure 7 as functions of the ratio 

Figures 7 (a) and 7(b) give the flutter boundaries in terms of the stiff- 
ness parameters bo^/V and boj^^V, respectively, and figure 7 (c) gives 

the frequency ratios a>ja>^ and atja^. In figures 7(a) and 7 (b) a rather 


sharp peak of required stiffness appears near a frequency ratio of unity 
when the structural -damping coefficient g is zero, as shown by the 
solid-line curve. Hie peak is greatly reduced if g is increased to 
only 0.03 (shown by the dashed line) and a further reduction is to be 
expected for higher values of g, but even large amounts of damping have 
relatively little effect on the required stiffness for frequency ratios 
less than about 0.6 or 0 . 7 . In actual structures, damping greater than 
the equivalent of g = 0.01 or 0.02 may be obtainable only by the addi- 
tion of damping devices. Figure 7(c) shows the frequency ratios 


and u>/ a >2 as functions of caija^. Again the solid-line curve is for 
g = 0 and the dashed-line curve for g = 0.03. For values of cu^ycog 

and 


less than 1.0 the flutter frequency et> falls between 


As 


<*V“2 


increases toward a value of 1 . 0 , 


a) is seen to approach both 


and and as a value of 1.0 is surpassed, tu drops below both co^_ 


o>l 

and 


0 ^ 2 . Structural damping has an almost negligible effect on the frequency. 


In the final section of appendix C the effect of using strip theory 
(multiplying the spanwise distribution of force and moment for a rigid 
translating wing by the flapping mode shape) to obtain lift and moment 
distribution for only flapping motion is examined. The resulting flutter 
boundaries are very close to those of figure 7 0 -ncl have, therefore, not 
been added to the figure. Thus, even though there are appreciable dif- 
ferences between the spanwise distributions of force and moment due to 
flapping motion as obtained by the two methods, virtually no difference 
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in the resulting boundaries was found for the particular parameters and 
modes of the present example. No generalization from these results is 
believed Justified, however, regarding the use of strip theory. 

In assessing the relative generality of the results of figure 7 , a 
fact to be kept in mind is that each of the two coupled modes of the 
control surface is held constant on the figure. For an actual structure, 
if the frequency ratio of the first two modes were changed by altering 
the flapping or pitching stiffness, or both, the coupled mode shapes them- 
selves would probably change and calculation of a new flutter boundary 
would be necessary. 


CONCLUDING REMARKS 


On the basis of linearized supersonic potential -flow theory, expres- 
sions have been developed for the velocity potential and for the section 
and total forces and moments on thin wings of arrowhead plan form with 
all sonic or supersonic edges oscillating harmonically in symmetric 
flapping and in roll. These force and moment expressions were derived 
by using the velocity potential expanded as a power series in terms of 
the frequency of oscillation and retaining terms to the third power of 
the frequency. The resulting total forces and moments are shown to be 
accurate over a range of frequencies wide enough to make them useful for 
many practical flutter applications. Multiplication of the section 
forces and moments of a rigid translating wing by the flapping and 
rolling mode shapes (called a "finite -wing strip theory") results in an 
overestimation of forces and moments in comparison with the results of 
the present analysis, which includes spanwise induced effects of varia- 
tion of section amplitudes. Use of the forces estimated by strip theory 
for the flapping motion in a single sample flutter solution, however, 
gave flutter boundaries virtually the same as the boundaries obtained 
with the forces given in the present analysis } however, no generaliza- 
tion is believed Justified regarding the use of strip theory. 

A method of flutter analysis using natural (coupled) vibration 
modes has been specifically adapted to an all -movable control surface 
configuration wherein all bending and twisting flexibility is concen- 
trated in the supporting shaft. With this method an example flutter 
analysis using two natural modes of vibration is carried out. A near 
equality of the first two natural frequencies was found to be very 
undesirable as regards the stiffness required to prevent flutter. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 8 , 1957* 



MCA IN 4189 


25 


APPENDIX A 

DEVELOPMENT OF VELOCITY POTENTIAL 


The characteristic coordinates used in reference 2 are also used 
in the present analysis, but in a nondimensional rather than a dimen- 
sional sense. (See fig. 3(b).) Thus, equation (4) can be written 



where 

u = |^(i - pn) 

"o = |^ x - 

v = ^-(i + PTl) 

v o = |^( x + Py) 

and w* is the downwash in transformed coordinates shown in figure 3( b )« 
As can be seen from this figure, the limits of integration in equa- 
tion (Al) are 
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where 

1 - pc 

' 1 + pc 


Symmetric Flapping Wing 

For the symmetric flapping mode of oscillation, the downwash 
w*(u,v,t) in equation (Al) is 


v*(u,v,t) = |^ f |v - u 


(A3) 


Upon substitution of equation (A3) into equation (Al), performance 
of the indicated integrations , and conversion to xy-coordinates , the 
result can be put in the form given by equation (6a). The coefficients 
of this equation can be expressed as follows: 


Af = a lf x + i2k (a 2 fX 2 + a,j f y 2 ^ + 4k 2 (a^fX^ + 

a^ f xy 2 j + i8k5^ag f x^ + ayfX 2 y 2 + a-g^y^ 

A f * = eqfy 2 + 4k 2 e 2f y^ 

B lf (x,y) = p lf (Cx + y) 2 + PgfP^yCCx + y) + 

i2k[j3jf (Cx + y) 5 + P 4f p 2 Cy(Cx + y) 2 ] + 
4k 2 [p 5f (Cx + y) 1 *- + P 5 f p 2 Cy(Cx + y)^| + ► 
i8k 5 j~P 7f (Cx + y) 5 + p gf p 2 Cy(Cx + y)^J 

B 2f (x,y) = B lf (x,-y) 


(A4) 

(A5) 


(A 6) 
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where 

a lf 

a 2f 

ajf 

a 4f 

®5f 

a 6f 

“7f 

a 8f 

e lf 

e 2f 

Plf 

P2f 


c£ 

a 


^|(2ct + 3) 


+ 3) 
3 a 2 


^SfL lM 2 (6a 2 + 22a + 15) + 2a 2 + 3al 
24pM L 



, [m 2 ( 24 o^ + 156a 2 + 235a +105)+ 24 a 5 + 75a 2 + 
36 op 6 cA L 


M^C 2 



+ 


283a 2 + 900a + 630) + 27 a^ + 270a 2 + 


= — Im 2 ( 4 a 5 - 4 a 2 + 65a + 105) - 6a 5 + 15a 2 + 45^1 

3603 2 a^ L J 
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03f - 


p 4f = 


P 5f = 


M^C 
2a^fa 

M 2 

2 p 2 0^ 

M 2 

I6p 2 a3/^ 


jM 2 (4cr + 5 ) + <yj 


Pfit - f Tp \f& + 5) + ^ 

12^Ca z Ja i - - 1 


3 7f 


Pflf = 


[m 2 (4c + 7) + 3^ 

48 P 2 a^ L J 

-^_[ M 2(2 0 + 5) + 3^ 


wherein 


a = p^ 2 - 1 


At the limiting condition of the sonic leading edge, regions II and 
III of figure 3 no longer exist and the integration of equation (Al) need 
be carried out only over region I. Thus, the coefficient A]_f of equa- 
tion (7a) is found to be 

A lf = T lf x + i2k 2 ^2fX 2 + Tjpy 2 ) + 4k 2 (r^f>x5 + T^xy 2 ) + 

i8k5 + T 7fX 2 y 2 + TQfy^) ( A 7) 


where 
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T 2f 

T 3f 

r kf 

T 5f 

T 6f 

T 7f 

T 8f 


-2M 2 

15P 4 


_ -SM g 


15P 

-M 2 

4203 ^ 

2 


_ M‘ 


8403 ^ 

9A538 

6 


_ -M 


3153 


6 


(13M 2 + 7) 
(13M 2 - 21) 
(5M 2 + 9) 


_ 

9^5 


(&4 2 - 9) 


Wing Oscillating in Roll 

For the wing oscillating in roll, the integrations of equation (Al) 
are carried out by using the following downvash in transformed coordi- 
nates (fig. 3(t>)} from equation (5b): 

v*(u,v,t) = ^ <p(v - u) (A8) 

Upon conversion to xy -coordinates, the coefficient of equation (6b) 

can be written 
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A<p = a lcp y + i2ka 2( pxy + 4k 2 + a^f^j + iSk^o^^y + (A9) 

where 


“KP = 

a 2cp “ 
a 3cp - 
a 4<p = 
a 5cp " 

a 6cp = 


-C 

T 


2M 2 C^ 

cr 2 

- M% 2. „ & 2 ( 28g + 45) + 9<yj 
24p 2 a3 L J 

— [m 2 (i 4 ct + 15) + 2a 2 + 3ffl 
24(B 2 c3 l 

|m 2 (19a 2 + 56a + 42) + 15 a 2 + 18a] 
36pV+ J 

^ 2 ( 22q + ^2) + 5a 2 + l8c] 

3 6p 2 cT L - J 


For the condition of the sonic leading edge, the coefficient A-^ 
of equation (7b) can be expressed as 

A 1(P ■ T lcpy + 12kT 2cp x y + ^ 2 ( T 3cp x2 y + T 4q>y 3 ) + i8fc3 ( T 5q^y + T6cpxy 5 ) (A10) 
where 


II 

_ -4m 2 

T 2 cp “ -5 

15P 5 

T 3 «p = -=^(5M 2 + 7) 
^ 315135 

TW, =■ ^5-(4m£ - 7 ) 

^ 515P 5 

T 5cp - =^Uv? + 9) 
9^5P 7 

T 6 cp ” 31 J 1 ™? - 9 ) 
9^5P 5 
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APPENDIX B 

FORCE AND MOMENT COEFFICIENTS 


Section Coefficients 

Symmetric flapping wing .- The coefficients 1^ f and M^_ f of 

equations (13) and (14) are obtained by substituting equations (6a) or 
(7a) , respectively, according to whether the leading edge is supersonic 
or sonic, into equations (ll) and (12). The general definitions uti- 
lized in reference 2 are as follows : 


V = £»"VvV*i a - eV 5 


■if e _ oe+nH-n 
r mn “ tL 


mi„_-i x i + 


(p^y^xv^Cx-i + y) n cos^ ~~ + 

0(Cxi + y) 


(-^y^x^C^ - y) n cos" 1 S u l } ~ y 

<hm - f Xl x-s/x 2 - fs 2 ^ to 

J y/o 


_ 2 e+nH-rH-l 


(p2Cy) e f Xl x“(Ox + y) n cos _1 2 +^jr 

Jy/C P(Cx + y) 


P x 1 

Jy/c 


(-0^Cy) e l * x^Cx - y)cos“ -L dx 


,-1 x - p^y 
0(Cx - y) 


(Bl) 


The following general definitions are also introduced: 


Jan - ~ ^ 

J x i + s|y| 

3^ - 2”*- n+1 f Xl - l EjM dx 

Jy/C Jx+P|y| J 


(B2) 
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With the general expressions of equations (Bl) and (B2) the coefficients 
1^ f and M ± f (i = 1,2) of equations (15) and (l4) for the condition 

of the supersonic leading edge can he written as follows: 

L 1 f ■ ^rif G 01 + a 2f E 02 + a 3f E 20 + e lf s 20 + Plf E 02° + ( B 2f H 0l" L + 

2rt 

P 5f F 03 ° + 3 4f F 02 1 + 4 2 (a 4f G 05 + a 5f G 21 + agfE^ + OffEgg + 
a 8f E 40 + e 2f s 40 + P5f H 04° + + p 7f F 05° + l 3 8f F 04 1 )| 

L 2,f =. 7^ ^-a-ifEoi " € lf J 20 “ 01f E O2° " f 3 2f F Ol’ L ) + k ( a 2f G 02 + 
a 3f G 20 " a 4 f E 03 “ a 5f S 21 " e 2 f J 40 + ^ 3 ^ 03 ° + P 4 f H 02 1 “ 

P5f F 04° “ P6f F 03 1 )| 


(B3) 


M l,f' = | a lf G 02 + a 2 f ( E 03 - G 02) + a 3f ( E 21 - G 2 o) + e lf S 21 + 

P If % 2° + P2f H ll 1 + Pjf ( E 13 ° _ H 03°) + P 4 f (^lS 1 ” H 02 1 ) + 
k2 ^4f G 04 + a 5f G 22 + a 6f ( E 05 " G 04) + °7f ( E 23 " G 22) + 
a 8f ( E 4l “ g 4o) + e 2f s 4l + l 3 5f H l4 0 + p 6f H 13 + 

P7f( F 15 ° " h 05°) + ^8f ( f 1 4 1 " E 04 1 )j| 

M 2 ,f’= ^ |^lf( G 01 - E 02) + e lf( S 20 - J 2l) + Plf^ 0 - F 12°) 
P 2f (Hqi 1 - F-jj 1 ! + k [agf G 03 + a 3f G 2L + % f (g Q5 - E o4 ) + 
a 5f ( G 2i - Egg) + egf (s 4o - J 4l) + P3f H 13° + P 4 f H 12 1 + 

P 5 f(Ho4° - F l4°) + - F 13 1 3} 




> (B4) 


+ 
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where 1 and M^j. are related by the expression 


M l,f + m 2,f = M l,f ' + m 2,t' ~ + iL 2,f) 


(B5) 


The coefficients $±f> and £j_f are defined in appendix A 

following. equation (A6). 

For the condition of the sonic leading edge the coefficients L ± f 
and reduce to 

Ll ,f = §flf G 01 + T 2f E 02 + T 3f E 20 + e lf S 20 + k 2 ( T 4f G 03 + T 5f G 21 
T 6f E 04 + T 7 f E 22 + T 8f E 40 + € 2f s 4o) 

k( _T lf E 01 " 6 lf J 2 o) + k ( T 2 f G 02 + T 3 f G 20 “ T 4 f E 03 “ 

51^21 " e 2f J 4o)] 


Lo f = — 

2it 


> (B 6) 


Ml >f ' " j T lf G 02 + T 2f( E 03 ” G 02) + T 3f( E 21 - G 2 o) + e lf S 21 
k2 Jj4f G o4 + T 5 f G 22 + T 6f ( e 05 - G oh-) + T 7 f ( e 23 - ^ 22 ) -» 
T 8f( E 4l " g 4o) + € 2f s 4lj| 

|j[jlf( G 01 - E 02) + e lf( S 20 - J 2l| + k [j 


Mo f » = •=! 

2jt 


r (B7) 


l T 2f G 03 + 


T 3f G 21 + T 4f ( G 03 “ E 0l) + T 5f ( g 21 “ e 22) + 
€ 2f ( s 4o “ J 4lj|j 


Equations (B6) and (B7) are obtainable from equations (B3) and (b 4) by 
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substituting for and dropping all terms involving the 

multiplier The coefficients ej_f and are defined following 

equations (A6) and (A 7 ), respectively. 

The main text of this paper includes a procedure for analyzing the 
flutter of wings or control surfaces of arrowhead plan form when the 
flutter involves vertical translation, pitching, and symmetric flapping. 
Equations (26) show that not only are section lift and pitching moment 
required, but also the product of the absolute value of the span coordi- 
nate and the section lift. This latter quantity for symaetric flapping 
is 

P|y'| = 2bPjyJ = -4pb 2 V 2 k 2 e iast f 0 (2|y|L lif + i2|y|L^ f j (B8) 

The quantities 2{y(L^ f and 2|yjL2 f are obtained by multiplying 

both sides of equations (B3) and (B6) by 2|y| . A comparison of the 
resulting equation with equations (Bl) and (B2) shows that, on the 
right-hand side of equations (B3) and (B 6), the index m is replaced 
by m + 1 in the quantities E^, C^, J^, and S mi1 , whereas 

is replaced by F mri e + ' 1 -/ P 2 C, and by g mri 6+1 y f p g C . To be precise, 

this substitution gives the desired quantities on the right (y positive) 
half-span, but since the lift P and the product P|y| are symmetric 
about y = 0 they became known also on the left half -span. The general 

definitions of F (nn e+1 and H mn e+1 are as follows: 

C 1 - 2 » 2 w p » 

= 2 e + 1+1 » n (|3 2 Cy) e+1 x 1 m (Cx 1 + y) n cos -1 - 

(-eSoyi^VtCx! - yJ'W 1 

San 6 * 1 - 

= 2 e+a+»n (p£c y je+l | 1 + y) n oos -1 dx - 

cl y 

(-pScy)^ 1 / 1 x m (Cx - yl^-coB- 1 x ~ ggCy , dx 

Jy/c 
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The general definitions of eq.ua tions (B9) (with the distinguishing bar 
under F and H) are thus seen to differ basically from those of equa- 
tions (Bl) only in that the second term within each square bracket is 
preceded by a minus rather than by a plus sign. Thus, on the right half- 
span, for the supersonic-leading-edge condition, 


2 |y[ L i,f = ^ | lx if G ii + a 2 f E i 2 + ^f^o + e if s jo + p^^ifSoa 1 + 

p2fSoi 2 + p 3f £ 03 1 + P 4f^02 2 ) + k 2 [ a 4f G 13 + W + 
a 6f E lk + “7^32 + a 8 f E 50 + £ 2 f s 50 + ^ + 

p 6fSo3 2 + ^fEc^ 1 + p 8f^o4 2 )]| 

2 |y| L 2 ,f = ^|^[- a if E n - e if J 30 + ^ ("Pii^oa 1 “ p 2 f^oi 2 ) + 

k [ a 2f G 12 + a 3f G 30 - a 4f E 13 * a 5f E 31 " e 2f J 50 + 

^'( p 3 f^ 03 1 + p 4f§02 ~ |3 5f^o4 1 " p 6f^03 3} 


V (BIO) 


For the condition of the sonic leading edge the quantities 2[y|L 1 ^ f 

are obtained from equations (BIO) by dropping all the terms that include 
a multiplier and replacing by T±f- The terms that include 

a multiplier remain unchanged. 


When the section coefficients are calculated, the chordwise inte- 
grals 0^, V e , ^ m e+1 j and Smn are needed. Most of the required 

integrals are given in appendix B of reference 2 or can be directly 
deduced from the spanwise integrals of equations (A23) and (A24) of ref- 
erence 5 . A few printing errors in reference 2 are noted in reference 5, 
and an errata has corrected equation (A23) of reference 5 as follows: 
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g 02 = f e 03 " %■ e 21 " ^P 1< ’ 1 4,4 


(Bll) 


Certain additional chordwise integrals needed for the flapping wing 
are as follows : 

Hop 0 - ^ *06° - ^[ 100 ^25 + 5 ° 2(30 + 11)El n 
64 (15a 2 + 70 a + 63)y^cosh- 1 ^ 

Hik 0 ' = -=±r Fq 6° + *15° “ -^[580^23 + C2(57a + 113)%! + 

■ -<2 w 5C ^ 120C 2 L 


30 C" 


64(57 a 2 + 170cr + 105)y^cosh _1 
H02 1 - jj - ^|c 2e 21 + l6(<1 + Sly^osh - 1 Sj 


> (B12) 


T, 1 _ 1_ V 1 P 2 v^ 

H 04 “ 5c F °5 20 


20^23 + 5 C 2 ^ a + 9 ) E 4 l + 


64 (3a 2 + 30a + 35)y 6 coslT 1 

,^ 5 -up ^ t(a,t3) > 


1 „ 2 P 

12 ” " ^2 "'04 ” 7^5 03 


Hno 1 = ~ F^. 1 ~ — F, 


v - ;th f <* 2 - M 20 ^ + 502(0 + 4)E u + 


Hl 5 pc 2 " U3 4p 2 c 3 


64(3ct 2 + 150 + lOjy^cosh -1 


Ohe chordwise integrals contained in equations (BIO) are 

defined by equations (B 9 ) and thus can be directly deduced from Bjjjn • 
For example, equation (A 24 l) of reference 5 gives Hqj , from which the 
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chordwise integral H q -^~ is obtained by merely removing the bars over 
the E's and the F. Then H is deduced to be 



_ 1 -pi 2 

" 4C ^04 


- 
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(B13) 


Similarly, from equation (A24c) of reference 5 ^ eciuce< ^ 4o 1:36 


So ? 1 - h V - hss2 fi W e ji + 2(50 - 5^'°^ £ 


X- 

Fy 


(B14) 


Consequently, a list of all the expressions for B am e+1 appearing in 
equations (BIO) is not needed herein. 

The remaining chordwise integrals needed in equations (B3) , (B4), 
and (BIO) are represented generically as follows: 


^nO = | Eml - 2i m+1 p2y IIH '^cosh -1 ^ 

^ 2 - P 2e dh-2,o) 

%SL " | ^3 - ^ W,! - ^V^cosh- 1 ^ 

*^n3 = 5 ®m4 " ^ hh - 2,2 ~ ^ ^^^^2+4,0 

SjnO = Jml - jj- Emo 

Sml = | J m2 - | Eml - 2m-lp2 ym f2 cosh -l ^ 


(B15) 


Other general expressions that can be useful for calculations are 
the reduction formulas 
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*mn “ p c ^ F m,ri+l 

w 

F m!> 2 - 

F 2 = 8 ^C 2 (f 0 

%i p \=m,iM-2 

Rolling wing.- The coefficients L^qj and Mj_ q, of equations (15) 

and (l4) are obtained by substituting equations (6b) or (7b) , according 
to whether the leading edge is supersonic or sonic, into equations (11) 
and (12). With the use of the general definitions of equations (Bl) and 
(B9) the coefficients may be written, for the condition of the supersonic 
leading edges, as follows: 

N 

L l,cp * 2“ a lq> G 10 + a 2cp E ll " P lf§02 “ P 2f5oi “ P 5f^05 " P 4f^02 + 

L. 

k2 ( a 5cp G 12 + a 4cp G 50 + a 5cp E 13 + a 6cp E 51 “ p 5fSo4° ~ ^6^0^ ~ 
p 7f^05 “ p 8f^04 ) * (BI7) 

^cp “ f“ |(- ct lcp E 10 + P ifEo2 0 + P 2f^l 1 ) + k ( a 2cp G ll " a 3q&-2 - 
a 4cp E 30 " p 3f5o3 " PltfSo2 + p 5f^04 + p 6f^03 ) 

M l,cp' = £ jv^ll + “^cp ( E 12 - ^l) - p lf^2° " P 2f&Ll 1 " 

p 3f (-13 “ §03 ) ” P 4f (-12 " 5q2 ) + k2 [^39 G 13 + 

a 4<p G 3l + a 5cp{^l.k ~ g 13) + a 6cp( E 32 - G 3l) “ £5^4° " 

' 

p 6f®L3 1 ” P 7f(— 15° " §05°) " P 8f(-l4 L " ^CA- 1 )" ' 

(Eq. (B18) concluded on next page) 



- cr M.l,n°) 

- ^-l.n 0 ) 

> (bi6) 

- + c2 V2,n°) 

- 2C ^l,n + l° + C ^2,n°) 
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M 2,<p' 


‘ a {^(“io - %) + - Eos 0 ) + 

!3 2f(-ll 1 " ^Ol 1 )] + k P‘2cp G 12 + {^12 ~ E 15) + 

a 4<p( G 30 " E 3l) ” P 3fSi3 0 “ ^f^ 1 + P 5f(-l4° ~ ^04°) + 

p 6f(43 L “ H 03 1 || 


(Bl8) 

Concluded 


where and <p are related by the expression 

M l,9 + ^2,cp = Mi,* 1 + - 2 x 0 ( 1 ^.^ + iLg^) (B19) 

For the condition of the sonic leading edge the coefficients L[_ 
and are obtainable from equations (BI 7 ) and (Bl8) by substituting 

r if for and dropping all terms involving a multiplier 


In the section contribution to the rolling moment of equation ( 16 ) 
the quantities 2 yL]_^ and 2 yLg jC p are obtained by multiplying both 

sides of equations (B17) and (Bl 8 ) by 2y. On the right side of these 
equations this procedure is equivalent to replacing by E^^n* 

V ^ <Wl,n, 4 n e ^ '■a** 1 /*® 0 ’ ^ ^ 1 W Wl /^' 

quantities are presented for completeness as follows: 


^ljCp = §j“ j a lcp G 20 + “^cp^gl ~ ■^(Plf H 02 1 + P2f E 01 2 + P3f E 03 L + 
04f F O2 2 ) + k2 |°'39 G 22 + a 4<p G 40 + a 5<p s 23 + “6(^41 “ 
-i-(p 5 fH 0 4 :1 ' + P6f H 03 2 + 07f F O5 1 + 08f F o4 2 )j| 

^ L 2,cp = |^|j a lcp E 20 + ^{Plf F 02 1 + 02f F Ol 2 j| + k [ a 2cp G 21 " 
a 3cp E 22 " a 4<p E 40 " , ^'(03f H O3 1 + 04f H O2 2 + P5f F o4 1 + 
P 6f F 03 2 )]} 


(B20) 
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For the condition of the sonic leading edge the quantities 2yL^^ are 

obtained from equations (B20) by dropping all the terms that include a 
multiplier 3^ and replacing aj_f by in the remaining terms . 


For calculation of 2yL^ C p and 2yLg ^ the following integrals, 

not included in equations (B12), (B15), nor in references 2 and 5 , are 
needed: 


H 05 1 " 35 f o 6 1 - ^^[^24 + C 2 (ka + 5 4)E 42 + 

(8c 2 + 1160 + l 83 )Egq] 

H 14 1 = ^2 F 06 1 - F 05 2 - ^[3C ^ 24 + C 2 (4a + 30)E 42 + 

(8c 2 + 68 a + 63 )E 6 q] > 

H 01 2 = 43 4 C 2 y 2 H 01 ° = ± F 02 2 - l6^Cja yVosh' 1 ^ 

H 02 2 = 4|3 4 C 2 y 2 H 02 0 = i. F 03 2 - | 3 4 C 2 ^ E 4 q 
Hq 5 2 = 43^2 y 2H 03 0 

=-^ f 0 4 2 '- £ + 64(30 + 5)y 6 cosh" 1 ^ 


(B21) 


If the section quantities of equations (Blj) and (Bl8) are desired, 
the quantities H^ 6 can be obtained from 1^® given in equations (B12) 

and (B21) herein and in reference 5 By use of the relations 


Smn 


0 _ 


2p 2 Cy 


•®mn 


H^ 1 = 2p 2 CyH flin 0 


(B22) 
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Thus, for example, in the first of equations (Bl8) can be deduced 

from Hq^ 1 in eq.ua tion (B21) to be 


TT 0 — 1 IT 0 

^05 -fe^06 


- 30^4 + c 2 (ka + 5^)e 32 + C 8 c 2 + 116 a + l83)E 5 q] (B 23 ) 


Total Coefficients 

The total lift force is obtained by substitution of equation (B 3 ) 
(or equation (B 6 ) for the sonic leading edge) into equation ( 17 ) . Total 
pitching moment is obtained by substitution of equation (B4) (or equa- 
tion (B7) for sonic leading edges) into equation (l 8 ). Root bending 
moment is obtained by substituting equation (BIO) into equation ( 19 ) . 
Total rolling moment is obtained by substituting equation (B20) into 
equation (20) . 


In order to obtain the spanwise integrals of E^, F, 


mn > 




*mn e . 


J mn J 


^mn) 


and H e , it is convenient first to reduce 


^mn y ^mn^ 


-mn ’ -mn 

in terms of 


Emn* 


■W COBh ' 


F e . 
-mn 


. -mn. — V in terms of F^, cosh" 1 and 

This reduction is obtained by use of equations (A23) and (A24) of refer- 
ence 5 and by means of the appropriate equations in the present appendix. 


Formulas for obtaining the spanwise integrals and 


F 6 
^mn 


are 


given in the appendix of reference 5* (Equation (Al4) of reference 5 


been corrected in an errata 

fie 

C- a iA. I 

should be 
IK 


The expression multiplying the quantity 


jEL 


2,n 
from 


(n + 1)(1 ± CD) 


The integrals 


F ® 
-mn 


are obtained 


F e 
x mn 


merely by changing the sign from plus to minus on the second 


major term (the term -which involves the double -primed quantities) in 
equations (AI 7 ) of reference 5- 


The only other integral needed Is 


J__ = 2 R . P . 

Tim 


vt 


2 mwl y m (l + Dy) n tanh - 1 1 + L B - ~ ftfr 

J 1 + (D + p)y 


d y 


(B2k) 


where R.P. means the real part of the succeeding quantity . By means 
of integration by parts 
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J m0 " 2 



2 m y m tanh“ 1 


1 + (D - P)y 
Ji + (d + Fly 


dy 


2* P 1 /^ ^dy 

m + 1 JO ^1 + 2Dy + (D 2 - B^y 2 



m 


(B25) 


where l4 jin is the I4 n of equation (A19) in reference 5« If the 

expression (l + Dy) n of equation (B24) is expanded, J r TTin can be 
obtained in terms of J m Q or I4 m . The needed recursion formulas are 
then 


J ml ~ 2J mO + ^Wl^O 


J m2 ■ 4j m0 + ^1,0 + ^m-2,0 


(B26) 


For the special case of the delta plan form (A^g =0), the total 

lift and pitching -moment coefficients due to flapping and the total 
rolling -moment coefficients due to roll reduce to concise forms, as 
follows i 


L 1 f = ~ - !^£(4M 2 + l)k 2 + + 12M 2 + 1)^ - 

' 3J3 5 45B 7 84O0 17 

(64M 6 + 24CM J+ + 12CM 2 + 5)k 6 

115,400p 15 

— M 6 C 2 (^4 + 2 CM 2 + 5)k5 
5240p i5 


> (B27) 
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* 

§c£_ _ ja fe? (]rti 2 + X ) k 2 + (^4 + 1SM 2 + 1 ^ k 4 _ 

150 5 10507 9450 11 


— ^ 6 ° 2 — (64m 6 + 24(34^ + 12QM 2 + 5)k 6 - 2xnl« r 
62570p 1 5 j 



^f.^! k + ^(1^2 + 3)k 5 . 

0 k 305 7209 

- M - c2 - (^4 lf + 20M 2 + 5)k5 - 2 xqLo r 
18OO0 1 5 *’ x 


> 


/ 


iL 1 . 4s2_ . S&W + 1)k 2 + _mV (8m 4 + lal 2 + 1)k 4 

150 5 3150' 189O0 11 

M^c5 — (61^6 + 240M^ + 12CM 2 + 5)k 6 

3118500 L 5 

k 

%,z - fs£ 1 . a&Z k + Mites (llM 2 + 3)k5 . 

^ 5P k 1505 2520^ 

— m6g3 - (8 m 4 + 20M 2 + 5)k5 
8100013 


(B28) 


(B29) 


The first two terms in. each of these total coefficients were obtained by 
spanwise integration of the section coefficients and also by the reverse- 
flow theorem as a check. The remaining terms containing higher powers of 
k were obtained by the reverse -flow theorem (using the velocity poten- 
tial for two-dimensional supersonic flow given by eq_. (23) of ref. 8) 
and were used in obtaining figure 5. 

From equations (B27) to (B29j_ the leading-edge sweep angle is seen 
to appear in such a manner that P/C^, M^/C 2 , and (Rolling foment )/c3 

are independent of the sweep angle of the leading edge for symmetrically 
flapping and rolling delta wings with supersonic leading edges^ This 
result corresponds to the finding in references 2 and 3 that'- P/C and 
BJC are independent of A for the supers onic -leading-edge delta "wing 

oscillating in vertical translation and pitch. Such simple dependence 
on A does not occur for the spanwise integrals of equation (19). 
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The total coefficients of equations (B27) to (B29) are plotted in 
figure 5 as a function of frequency for a Mach number of 1.6j a dis- 
cussion of figure 5 is given in the section entitled "Accuracy as a 
function of reduced frequency." 
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APPENDIX C 

NUMERICAL EXAMPLE OF A FLUTTER CALCULATION 


Calculation by Present Method 

A numerical example is now given of the method of flutter analysis 
outlined in the main text. The configuration chosen (see fig. 6) iB 
that of a half -span all -movable control surface of arrowhead plan form, 
with A = 45° and A^g = -15°, mounted on a shaft and with the wind- 

tunnel wall presumed to act as a reflecting plane. The control- surface 
pitch axis is located at Xg = 0.55^ (40.5 percent mean aerodynamic 

chord) . The location of first-mode and second-mode node lines are given 
in terms of the nondimensional coordinates (x,y) and these node lines 
are straight because all flexibility is considered to be concentrated 
in the supporting shaft. The components of deflection, normalized to 
a tip vertical deflection downward equal to the root chord 2b, are 
as follows: 

h, 

= 0.0600 
2b 

9 X = 1.071 
f-L = 0.855 


^2 _ 
2b 

0 2 = 
fo - 


0.232 
- 6.02 
2.77 


(Cl) 


The Integrals over the half-span model of the necessary mass properties 
are as follows: 


Jm dy = 0.00253 slug 
S a dy = 0.0000211 slug-ft 

dy = 0.000232 slug-ft | 

J\ a dy = 0.0000227 slug-ft 2 
J 1 If dy = 0.0000433 slug-ft 2 
/ p af dy = O.OOOOO855 slug-ft 2 J 


(C2) 
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Insertion of these quantities in equation (22.) gives the gener al ! zed 
mass terms 


Aj_ = 0.0000880 slug -ft 2 
Ag = 0.001014 s lug -ft 2 


(C3) 


and with b = 0.2375 and p = O.OOO66, 

a L = 12.04 


a 2 == 254.4 


For a Mach number of 1.6 values for the aerodynamic half -span 
integrals of equation (26) are as follows: 


r ± = 0.10695 - 0.11813k 2 + 1(0.30190 J - 0.12790k) 
r 2 = 0.041546 - 0.033206k 2 + i (0.093984 J - 0. 040931k) 

I3 = O.30190 - 0.14464 + i (-0.07825 ^ + 0.16227k) 

k 

r 4 = 0.093984 - 0.053708 + i( -0.034456 4 + 0.05l46lk) 

k 

r 5 = 0.048896 - 0.035692k 2 + i(o. 16435 4 - 0. 046946k) t- 

r 6 = 0.016912 - 0.0045837k 2 + i (0.064980 4 - 0.013625k) 

r ? - 0.01884 - 0.o4704k 2 + i (0.01525 i - 0.04004k) 
r Q = 0.01525 — - 0.029437 + 1(0. 04448 i + 0.05277k) 

r 9 = 0.014556 - 0.023179k 2 + 1(0.02824 4 - 0. 018644k) 


(C4) 
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The aerodynamic integrals C ^n of equation ( 25 ) for the coupled 
modes are specified by equation (Cl) are then 

C^L = - 0.22841 + 0 . 09515 ^ + i(- 0 .l 4 ll 8 | - 0 . 11845 k) 


(C5) 


c 12 = 1.2839 - 1.03499 + 1^-0.66961 ^ + 1.1822k) 

= -0.60934 ^ + 0.15620 + i (0.23419 ^ - 0.053134k) 
c 22 = 3.4259 \ - 1.2857 + i(- 3 . 98 l 6 i + 0.23698k) 


The values of a^ and Cm-n listed above are applied in the flutter 
determinant (eq. (28)), and the determinant is solved over a range of 
values of k for the two real unknowns cd^cd and o^Jcn. Of course, 

g L and g 2 can be varied independently, but in the present example the 
simplification gj_ = g 2 = g is made. The flutter results, Bhown in 
figure 7 as functions of the frequency ratio a^/c for two values of 

the structural damping, coefficient g, are discussed in the body of the 
report in the section entitled "Discussion." 

Effects of Using an Aerodynamic Strip Theory 
for Flapping Motion 

In the discussion of figure 4 in the main text, the differences in 
the section lift and moment coefficients as obtained from the expressions 
of the present analysis and as obtained alternatively by a "finite-wing 
strip theory" were compared and found to range frcm moderate to appre- 
ciable. In order to determine what effect the use of such strip-theory 
coefficients has on the present flutter example, substitutions can be 
made in equation (26) as follows: 


In and 1 ^: 


Use 2 |yj^Lj_ + iL 2 ) for ^L].,f + iL 2 ,f) 


(c6a) 
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In r 9 : 


Use 


y|(M 1 + iM 2 ) for (M ljf + lH 2)f ) 


( C6b ) 


As a result of these substitutions, 

r 5 = 0.083092 + 1(0.18797 i - 0. 08l862k) 
Tg = 0.039631 + i (o.082213 £ - 0.035867k) 
r 9 = 0.018949 + i (0.028100 i - 0.025442k) 


(C7) 


The other values are as listed in equation (c4). A brief discus- 
sion of the effect of using an aerodynamic strip theory to account for 
flapping motion in the present example is also Included in the section 
entitled "Discussion" in the main text. 
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(a) Section lift coefficients for symmetric 
flapping mode. 


(t>) Section pitching-mcment coefficients for 
symmetric flapping mode. 


Figure 4.- Spanwise distribution of section lift and pitching -moment coefficients for symmetric 
flapping mode and of lift coefficients for rolling mode for a V?° delta wing (C = 1, D =■ 0) 
with Xq - 0.5; M - 1.6, and k - O.p. Shown is the effect of using one and two terms of 


the frequency expansion for the flapping and rolling wing and of using two frequency terms 
from a finite-wing strip theory. 
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(c) Section lift coefficients for rolling mode. 
Figure Concluded. 


and 




d- 0) 













Figure 6.- Schematic diagram of arrovhead-plan-foim all-movable control surface analyzed in 

appendix C . 
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(a) Stiffness parameter, 



boj, 

(b) Stiffness parameter, 


Figure 7.- Flutter boundaries in terms of tbe stiffness parameters 
bcnp/v and bo^V as functions of tbe natural-mode frequency- 

ratio 0^/002 for two values of the structural damping coeffi- 
cient g an d the associated flutter frequency ratios oyb^ 
and cd/o) 2 at M = 1.6 for the configuration sho>m in figure 6 
and analyzed in appendix C. 
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(c) Frequency ratios 



and 



Figure 7*- Concluded. 
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